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The free space dose rates from the various components of the 
galactic cosmic rays can be determined from the particle energy 
spectra and energy loss data. Depth-dose profiles in material 
shields or in body tissue are complicated by the nuclear cascade 
set up in the material by the incident primaries. While consid- 
erable work has been done on the proton induced particle cascade, 
the heavy primary secondary dose from the galactic cosmic rays, 
which also fragment and cascade, has not been estimated. 

In this paper we present the results of a calculation which 
estimates the heavy primary secondary doses from cosmic ray inter- 
action data. The incident galactic cosmic ray heavy primary spec- 
trum is represented as the sum of helium, nitrogen (M group), mag- 
nesium (light heavy), and iron (very heavy) components. The inci- 
dent iron nuclei are allowed to fragment into lesser Z secondaries, 
which are assumed to travel in the same direction and start with 
the same energy per nucleon as the interacting primary. The total 
emergent particle energy spectra and dose are then presented for 
the galactic heavy primary spectrum incident on aluminum and tis- 
sue slabs. The importance of the fragmentation parameters assumed 
is also evaluated. 

The total dose from the heavy primaries and their secondaries 
is found to be reduced by only a factor of two in 20 g/cm^ of 
shielding . 

*Work supported by NASw-1963 


We present here a method for the cal- 
culation of the secondary particles pro- 
duced by heavy galactic primaries. The 
galactic cosmic ray spectrum consists of 
protons, helium and higher Z elements, all 
with a similar energy spectrum, but of re- 
duced intensity with increasing Z. As 
these particles pass through matter, both 
nuclear interaction and ionization reduce 
the intensity and energy of the incident 
spectrum. While the electromagnetic inter- 
action for the particle energies of inter- 
est are well understood, and their effects 
readily calculated, the nuclear interac- 
tions lead to complex cascades of secon- 
dary particles. While a great deal is 
known about the cascade products of pro- 
ton Induced reactions, the cascades in- 
duced by the higher Z elements of the gal- 
actic cosmic rays have received less atten- 
tion . 

The fraction of the yearly free space 
dose resulting from the various components 
of the galactic cosmic rays are shown in 
Table I, as determined from the particle 
energy spectra given in Reference 1. 

In Figure 1, we show the depth-dose 
profiles for the helium and higher Z par- 
ticle groups under two assumptions: 1) that 
only electromagnetic interactions attenuate 
the flux, and 2) that the particles suffer- 
ing nuclear interaction are removed from 
the beam (uncollided particle dose only). 
The wide range between these estimates In- 
dicates a need for further analysis. 


TABLE 

Component 

I 

Yearly Dose 

H 

(z- 1 ) 

4.6 

He 

(Z = 2) 

3.5 

M 

(6<Z<9) 

1 . 9 

LH 

(10<Z<14) 

1.3 

VH 

( 26 <Z <28 ) 

1.3 


TOTAL: 12.6 rads 


year 
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Figure 1: HEAVY PARTICLE GALACTIC DOSES 


X 


FORMULATION OF THE CALCULATION 

Let ^ (E , x) be the differential num- 
ber flux of ions of type i at a point x in 
the absorber. Then it is easy to show that 
the relation between the incident flux and 
the flux at a point in the absorber is 
given by 

* i <E’,x)=* i <E,0)f£f4 ) 


neglecting nuclear interactions. 

If we include nuc.lear interactions, 
and change variables to energy per nucleon, 
T , we have 


Cl) 


S (T) 
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,x) = <^ i ( T,0 )x 

^CDdt 

exp ^£ s (T) 


where S(T) is the stopping power of the ion 
of type i and energy per nucleon T. If M 4 
is energy dependent, then we have simply 
exp [“M-jx] for an attenuation factor. The 
uncollided flux can then be determined at 
points of interest in the absorber. 

Now ions of type i produce secondary 
ions of type j at points in the absorber 
at the following rate 


< 2) ^L(t 


x) = ^ i (T,x)^ a i (T) p i j( T, T') 


where 

(T , x) is the energy spectrum of ions 
of type i 

a^CT) is the interaction cross sec- 
tion for ions of type i and 
energy T 

P (T,T*) is the fragmentation parameter, 
J giving the probability for an 

interaction of type i ions at 
energy T to produce a secon- 
dary ion of type j and energy 
T ’ 

N q - Avogadro’s number 

M = the atomic weight of the ab- 
sorber 


d<£ 

-j^-(T f ,x)= the rate at which particles 
of type j and energy T' are 
being produced in the absorber 
at point x 


Particles of type j are slowed and 
attenuated as they pass through the re- 
mainder of the absorber and they emerge 
with an energy related to their initial 
energy T, and remaining path length in 
the absorber, X-x , of 
RJ (T) -Rj (T f ) +X-x 

and with an intensity reduction of 
f(777 exp[- Uj (X-x>] 


The emerging flux of type j particles pro- 
duced by type i primaries is then 


(3) (T f.X)"! 2 - J<p ± (T,x) 
0 

V T)P ij (T) tiff) x 


exp [-u j (X-x) ] dx 

Now we apply these basic relations 
to the problem of determining the final 
penetrating energy spectrum of the gal- 
actic heavy primaries and induced secon- 
daries in an absorber. First, we con- 
sider a maximum Z of 26, as only a small 
number of higher Z particles are present 
in the galactic cosmic ray flux. The 
penetrating flux of 2 6 Fe P art l c l es at an Y 
point in the absorber can be determined 
by expression (1). Once ( T »x) is 

known, then the rate at which lower ato- 
mic number particles are being produced 
at points in the absorber can be deter- 
mined bv expression (2). 


The total production rate of particles 
of type j is given by the following sum 

26 

ffi( T ,x)= X IrV^ijd) 

i«j+l 

where <p,( T,x), the total flux of particles 
of type i, depends on both the incident 
flux and the accumulated secondary flux. 

To determine the penetrating flux of 
particles, a computer program was devel- 
oped to evaluate the expressions discussed 
numerically. This program evaluates the 
production rates of the secondaries at a 
set of thickness points in the absorber for 
a suitable energy grid. Then by sequential 
transporting the highest Z element through 
the absorber and determining the production 
rates for all lower Z elements, we obtain 
the final penetrating number energy spec- 
trum . 


BASIC ASSUMPTIONS USED IN THE CALCULATION 

The preceding section has developed 
the formations used in the calculation of 
the secondary particles produced by the 
heavy galactic primaries. Now we discuss 
the basic assumptions used in this analy- 
sis, and some of the resulting limitations 
in the present calculations. 

First we employ the straight ahead, 
or one dimensional approximation. This 
can be justified by use of cosmic ray in- 
teraction data as discussed in Reference 2 
which shows that the heavy fragments emerg- 
ing from the interaction site are strongly 
peaked in the forward direction. It is 
well known that the primary heavy parti- 
cles suffer little angular deflection while 
passing through matter. 

Second, we use the energy independent 
overlap cross section model used by Cleg- 
horn, et al. (Reference 3) to fit his cos- 
mic ray emulsion results. It was found 
that cross sections were essentially inde- 
pendent of energy in the energy range 100 
MeV/nucleon to 30 GeV/nucleon. At energies 
below 100 MeV/nucleon, the higher Z pri- 
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maries have a small range, but their cross 
section will probably be underestimated 
in this energy region. 

The fragmentation function, P^.(T,T f ) 
represents the least well known quantity 
required in the calculation. Since the 
energy dependence of the secondary frag- 
ments has not yet been described, we make 
the simple assumption that T=T 1 , that is 
the fragments have the same velocity as 
the primary ion. This assumption is in 
agreement with an interaction model in 
which the heavy primary is stripped of 
some fraction of its mass, and proceeds 
on with unchanged velocity. Two sets of 
fragmentation parameters were then used 
to explore the effects of various assump- 
tions on the final results. In the first 
approximation, we have used a set of frag- 
mentation parameters which assume one 1 H, 
one ^He , and an equal probability of higher 
Z fragments normalized to conserve Z in 
the incident heavy primary. This approxi- 
mation, when coupled to the equal velocity 
approximation, leads to approximate con- 
servation of energy, neglecting nuclear 
binding energies and pion formation. The 
second set of fragmentation parameters 
are described in the paper by Curtis in 
these Symposium proceedings. 

Both sets of fragmentation parameters 
are assumed to be independent of the pri- 
mary particle energy, a conclusion that 
is consistent with the results of Cleg- 
ho rne . 

Finally, the calculations have neg- 
lected the interactions of primary and 
secondary ions with hydrogen. While the 
results of Bertini’s internuclear cascade 
calculation can be applied to this prob- 
lem by a proper transformation of rest 
frames this has not yet been incorporated 
in the calculation. The results given 
for tissue thus neglect the fragmentation 
induced by hydrogen. 

RESULTS 

Typical results of this c alcula t ional 
method have been developed with the assump- 
tions discussed. First, using the VH 
group heavy particle spectrum presented 
in Reference 1, depth dose profiles for 
a unidirectional beam of particles inci- 
dent normally on a slab of water and alu- 
minum are shown in Figures 2 & 3. The 
total tissue dose, the 2 6 Fe ( VH ) dose, 
and selected lower Z secondaries are 
shown. For the two sets of fragmen- 

tation parameters are represented by (1) 
solid lines and (2) dotted lines. We see 
the expected rise and more gradual fall 
of the secondary dose components. 


• Set I Fragmentation 
Parameter* 


-Set II Fragmentation 
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Figure 2: HEAVY GALACTIC PRIMARY VH GROUP DOSE 
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Figure 3: HEAVY GALACTIC PRIMARY VH GROUP 



With the complete GCR spectra, as 
represented by VH, LH, M and Helium com- 
ponents, the results are as shown in Fig- 
ure 4 . 



Meyer (Reference 4) has reviewed the 
galactic cosmic ray data available up to 
1969, and has presented an estimate of the 
intensity of all the heavy primary compo- 
nents below 26 ^ e * From this compilation 
we have calculated the total depth dose 
profiles in water for incident primaries 
2 He to fi Fe , and presented the results 
In Figure 5. 



Figure 5: HEAVY GALACTIC PRIMARY TOTAL DOSE 

Separate Spectra, Excluding Proton Dose 


CONCLUSIONS 

The results presented indicate the 
general nature of the total depth-dose pro 
files can expect in tissue and alumi- 

num absorbers. Uncertainties in the frag- 
mentation parameters lead to significant 
changes in the total dose, as well as the 
dose distribution among the secondary com- 
ponents. Further work in defining the 
fragmentation parameters and incorporation 
of the hydrogen interaction results for 
tissue would improve these dose estimates. 

While only dose results are presented 
in this paper, the penetrating particle 
energy spectra are also calculated, and 
can be used to examine the LET spectra for 
further study of the biological implica- 
tions of the galactic cosmic rays. 
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